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Abstract 

Parametric correlations of energy spectra of quantum chaotic systems are 
presented in the orthogonal-unitary and symplectic-unitary crossover region. 
The spectra are allowed to disperse as a function of two external perturba- 
tions: one of which preserves time-reversal symmetry, while the other violates 
it. Exact analytical expressions for the parametric two-point autocorrelation 
function of the density of states are derived in the crossover region by means 
of the supermatrix method. For the orthogonal-unitary crossover, the velocity 

distributions is determined and shown to deviate from Gaussian. 
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Spectral properties of quantum chaotic systems are well described by the Wigner-Dyson 
statistics of random matrix theory |L|. Originally introduced to study complex nuclei, ran- 
dom matrix theory (RMT) was found to be relevant to the spectral statistics of disordered 
metallic particles (quantum dots) @-[§], as well as quantum chaotic systems in general (e.g., 
chaotic billiards) The most striking feature is universality: when energies scaled by the 
mean level spacing, their distribution depends only on the symmetry of the Hamiltonian, 
irrespective of microscopic details. Typically systems belong to one of three universality 
classes: orthogonal (spinless) and symplectic (with spin-orbit interaction) for those which 
are invariant under time-reversal invariant systems, and unitary for those which are not. 
Starting from the Schrodinger equation with a random potential, Efetov has developed a 
field theoretic description of disordered metals, based on a supermatrix nonlinear cr-model [|J. 
He demonstrated that, up to the Thouless energy, the density of states correlators for the 
three universality classes are identical to those derived from RMT. 

Subsequently Pandey and Mehta extended RMT to examine crossover behavior where 
the symmetry of the system gradually changes from one pure symmetry to another ||. 
They obtained an analytical expression for the density of states correlator, i.e., the two- 
point cluster function, and related quantities. (A result recently rederived recently by the 
supermatrix method 0.) The resulting description of the orthogonal- unitary crossover was 
numerically justified for a disordered metallic ring with Aharonov-Bohm flux ||. 

A second approach to quantum chaos is to examine correlations of a Hamiltonian H(X) 
which depend on an external parameter X. Beginning with the work of Pechukas and 
Yukawa |10| |, problems of this kind have attracted great interest most recently |11| , 12| . By 



making use of an appropriate rescaling, parametric correlations of the density of states as 



well as response functions were also shown to be universal 

In this Letter, we will examine the universal parametric correlations in the orthogonal- 
unitary crossover region. For completeness, we mention (but do not discuss) the result 
for symplectic- unitary crossover. We are motivated by problems where the Hamiltonian 
H(X ,X U ) depends on two external parameters X Q and X u , where the former preserves 



the time-reversal symmetry whereas the second violates it. Although the spectral statis- 
tics of H(X Q ,X U = 0) belong to the orthogonal ensemble, gradually driven to unitary by 
introducing T-invariance breaking parameter, X u . 

Parametric correlations in the orthogonal-unitary crossover are relevant to many prob- 
lems. In a recent experiment flI5|, the correlator of differential conductance was measured 



as a function of magnetic field, X u , in a heavily doped quantum dot. These measurements 
were directly related to density correlators in the crossover region, considered below. At the 
same time, one can imagine gate voltage serving as second parameter X Q , which perturbs the 
system but conserves T-invariance. A second example could involve an irregular ballistic 
cavity, or billiard in which a potential, X , changes the shapes of the boundary, while a 
magnetic field acts T-breaking perturbation. 

Our aim here is to determine the parametric autocorrelator of density of states 

K(Q, X , X u - X) = (u{E, X - f , X - 4f ) v{E + n, x + fs x + fo) , (i) 

where v{E, X Q , X u ) = J2 n ~~ E n (X a , X u )) is the density of states and E n (X Q , X u ) are 
eigenenergies of the Hamiltonian H(X , X u ). The statistical average denoted by (• ■ ■) can 
be performed over certain interval of the energy and/or external parameters. K depends on 
the average of the unitary parameter X, as well as on differences of parameters X 0>u and 
energies, Q. X serves as a crossover parameter from orthogonal to unitary symmetry. 

Following Ref. p3| , |14|| , we switch from {Q,X ,X U ,X} to dimensionless variables 
{u, x Q , x u , x} (conventionally written as lower case), 

u = n/A] e n = EJA (2) 



x = v / C o (0;X = 0)X o ; x u = yjC u (0; X -»■ oo) X u ; x = yfC u (0; X -»■ oo) X, (3) 

where A = is the mean level spacing of the spectrum. We note that rescaling for X Q 

is defined in the orthogonal limit, and for X u and X at the unitary limit, which follows the 
conventions of both. We show that a dimensionless correlator defined by 



k(u, x , x u ; x) = A 2 K(Q, X , X u ; X) - I, 



(5) 



is a universal function of all dimensionless variables lu, x 0jU , and x. 

Consider the N x N random matrix H(X D ,X U ) = Hq + X Q H S + X u H a , where Hq and 
H s are real symmetric matrices and H a is an antihermitian matrix |J, the transition from 
orthogonal to unitary symmetry for K(VL, X , X u ; X) occurs at X of the order of 1/y/N or 
equivalently x ~ 1. Thus all non-universal effects such as X-dependence of A are negligible 
in the limit N — > oo. As a result, crossover behavior of k(u, x Q , x u \ x) is universal. 

To evaluate the k(u>, x Q , x u ; x), we use the zero dimensional supermatrix nonlinear a- 
model, which can be derived from the disordered metallic grain with large dimensionless 
conductance. This model serves as an underlying universal model for describing spectral 
statistics of quantum chaotic systems HT3| , |14| even in the orthogonal- unitary crossover re- 
gion j7|. Following the notation of Ref. |J, k(u be expressed through the 
nonlinear a-model as 



k(u, x u , x - x) = -Re J VQ e - F ^(Q - A)H(Q - A) 
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where Q is a supermatrix. The extension to parametric correlation functions is straight- 
forwardly done | 16|| . Performing the definite integral over the supermatrix, we obtain for 
orthogonal-unitary crossover: 



where 



k ou (co, x OJ x u ; x) = Re j dX J dXi J d\ 2 e W + (uj,x ,x u ;x), 
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(10) 



and we define a = tc 2 (x 2 — rr^/4)|l — A 2 |. For completeness, we include the analogous 
expression for symplectic-unitary crossover: 



roc rl rl 

k su (uj,x ,x u ;x) = Re / dX dXx / dX 2 e~ F W-(uj, x a , x u ; x), 

Ji J-i Jo 



(11) 



where x D plays a role as any T-invariance perturbation parameter in this case. Hereafter we 
restrict ourselves to the orthogonal- unitary crossover case for brevity. 

We first check various limiting behaviors of Eq. (||). The orthogonal limit straightfor- 
wardly gives the same result of Ref. when we set x, x u — > 0. In the unitary limit (x — > oo, 
and x Q — ■> 0), we can recognize the leading contribution coming from the term proportional 
to x 2 exp (-F + a). After A2 integration, the expression is reduced to the same expression 
for the pure unitary case in Ref. |13[]. Another interesting limit is the large uo asymptotics, 
which can be obtained by expanding the integrand of Eq. (f|) around A = Ai = A 2 = 1, and 
replacing the integral region over A by (—00, 1). After straightforward but lengthy calcula- 
tion, we can show that the large a;-asymptotics acquires the diffusive form, which can be 
interpreted as the Diffuson and Cooperon modes in disordered systems: 
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We can use k(u, x u , x Q ; x) to obtain the velocity (i.e., single level current de n /dx 0)U ) dis- 
tribution functions. Here we can define the two kinds of rescaled velocities both in the or- 
thogonal and unitary directions by de n /dx 0tU . Their probability distribution can be obtained 
from the formula f 0)U (v; x) = (6 (v - dej dx 0jU )) = lim^^^^o x ,u K u = vx o,u, x a , x u ; x) 
Substituting our exact result Eq. (|§|), we obtain 
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~ 2w2si2 Jq 1 e 2n2s:2x2 dX. The function A(x) behaves like 1 - 47r 2 x 2 /3 for small 



V2tt 

where A(x) = e 

x, but vanishes exponentially for large x. Their limiting behaviors are easy to understand. 
In the unitary limit (x —>■ 00), the velocity distribution has a Gaussian in both directions, 



(13a) 
(13b) 



ft - \ ex P h w f, - \ exp -w 2 /2 

f (v; x -»• oo) = = — ; f u {v;x^oo) = = . (14) 
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This means that there is no qualitative difference between the unitary and the orthogonal 
perturbations, since T-invariance is already fully broken in the unitary limit. On the other 
hand, in the orthogonal limit, the two velocity distributions behave quite differently. 

f (v; x^0) = 6XP [ ^! /2] ; f u (v;x^0) = 6(v). (15) 

V Z7T 

f u {v;x — > 0) collapses to a 5-function, since e(x Q ,x u ) — e(x o ,0) oc x 2 for x u — > 0. However 
there is nothing special in the orthogonal direction, so that its distribution is again Gaussian 
with twice as wide variance as the unitary limit. The velocity distribution interpolates 
smoothly between these two limiting cases. 

Since the average velocities vanish, the velocity distributions in the crossover region can 
be characterized by their variances: c OjU (0;x) = ((de n /dx 0tU ) 2 ). According to Eq. (11), they 
are 

. 1 + A(x) (2ir 2 x 2 + l)A(x) - 1 . 9 9 , , 9 9 , , 
Co (0; x)= V ; + - Ei(-2vr 2 x 2 ) exp(2vr 2 x 2 ), (16a) 

c u (0; x) = 1 - [l + 2vr 2 x 2 exp(27r 2 x 2 ) Ei(-2vr 2 x 2 )] A(x), (16b) 

where Ei(— z) = f£° dte~*/t is the exponential integral function. Hence in the limiting case 
for i < 1, we find 



c o (0; x -> +0)^ 1 + ln(2vr 2 x 2 ), (17a) 
c u (0; x -> +0)= -2tt 2 x 2 ln(27r 2 x 2 ). (17b) 

These logarithmic dependencies can be understood by a 2x2 random matrix model fTSfl, since 



the nearest neighboring pairs of energy levels make the dominant contribution for small x. 



Around the orthogonal limit, eigenenergies can be considered as e(x) ~ ye 2 (0) +x 2 , where 
e is assumed to obey the Gaussian orthogonal ensemble. Therefore c u (0; x) can be evaluated 

as 

c u (0;x) w l°° de (de/dxf Pgoe(^) ~ l_dex 2 /e « -x 2 ln(x). (18) 



This demonstrates that the logarithmic dependence is specific to the orthogonal-unitary 
crossover. By contrast, c n (0; x) ~ x 5 in the symplectic-unitary crossover because Pgse(^) ~ 
£ 4 for small e. 

In Fig. [L], we compare our analytical results of c D)U with numerical results, obtained from 
the random matrix H = H + X H S + X u H a , where H , H s and H a were defined previously. 
((dE n /dX 0tU ) 2 ) was determined numerically and rescaled according to Eqs. (|HI|)- As is seen 
in Fig. p], agreement with Eqs. (|16a| ,b) is good, particularly for small x. 

So far, we have used the rescaled parameters which are determined in the unitary or 
orthogonal limit. However, in experimental situations, only the results in the crossover 
region may be obtained. For such a case, we propose here the local normalization scheme. 
While x and x u are formerly defined by Eqs. @, what is actually observed at a fixed X is 

To express all universal properties in terms of y and y u instead of x and x u , we can make 
use of the relation between y and x {y , u an d x 0)U ) and write 

x 2 c u {x) = f- x 2 0)U c , u (x) = y 2 ou . (20) 

By calculating the function x = f(y) and x QjU = y ,u/c ,u(f(y)), we can express the universal 
dependence in terms of {y ,y u ,y} rather than {x a ,x u ,x}, which enables us to rescale the 
external parameters through a local point in the crossover region. 

In conclusion, we have presented an analytical expression for the universal parametric cor- 
relations of density of states both in the orthogonal-unitary and symplectic-unitary crossover 
region of quantum chaotic systems. For the former, the velocity distribution functions both 
in the unitary and the orthogonal direction were presented analytically. In particular, the 
velocity distribution in the unitary direction shows continuous transition between a Gaus- 
sian and a 5-function. Their variances both in the orthogonal and unitary directions are 
shown to have the logarithmic dependence around the orthogonal limit. In addition, we 
have also proposed a local rescaling scheme of data in the crossover region. 
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FIG. 1. Comparison between numerical and analytical calculations of c OiU (0;x). Numerical 
calculations of c DiU were obtained from 25 x 25(0), 50 x 50(A) and 100 x 100(v) random matrices 
H(X Q , X u ) = Hq + X Q H S + X u H a described in the text. The dot-dashed and the dashed lines are 
analytical results for c u (0;x) and c o (0;x) given by Eqs. ( 164 b). 
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